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We present a search for the decays B° — > e^e , B" — > fj,'^ fi and B° — > e*/i^ using data collected 
with the BABAR detector at the PEP-II e'^e~ collider at SLAC. Using a dataset corresponding 
to 384 X 10*^ BB pairs, we do not find evidence of any of the three decay modes. We obtain 
upper limit on the branching fractions, at 90% confidence level, of B{B'^ e^e~) < 11.3 x 10"*, 
B{B° ^ < 5.2 X 10"*, and 23(5° e^y."^) < 9.2 x 10"^ 

PACS numbers: 13.20.He,14.40.Nd 



The standard model (SM) of particle physics does not 
allow flavor changing neutral currents at tree-level, and 
decays of this kind are predicted to have very small 
branching fractions. This makes rare decays particularly 
interesting for the detection of possible new physics (NP) 
beyond the SM, such as supersymmetry 1] (SUSY): loop 
contributions from heavy partners of the SM particles 
predicted in these models might induce, for certain de- 
cay modes, branching fractions significantly larger than 
the values predicted by the SM. 

The leptonic decays 5° I'^l'^ (where stands 
for e+e~, or e^/i^; charge conjugation is implied 

throughout) are particularly interesting among rare de- 
cays, since a prediction of the decay rate in the context 
of the SM can be obtained with a relatively small er- 
ror, due to the limited impact of long-distance hadronic 
corrections In the SM, B'^ decays proceed 

through diagrams such as those shown in Fig. [T] These 
contributions are highly suppressed since they involve a 
b —f d transition and require an internal quark annihi- 
lation within the B meson. The decays are also helicity 
suppressed by factors of (mi/mB)'^ , where mg is the mass 
of the lepton and tub the mass of the B meson. 

In addition, 5° decays to leptons of two different fla- 
vors violate lepton flavor conservation, so they are for- 
bidden in the SM. This feature provides a handle to dis- 
criminate among different NP models 4]. 

The B'^ decays are sensitive to NP also in 

a large set of models with Minimal Flavor Violation [5] 
(MFV), in which the NP Lagrangian is ffavor blind at 
the typical mass scale of new heavy states, with reduced 
effects on flavor physics at the B mass scale Q. In the 
context of MFV models, NP corrections to B° 
are characterized by interesting correlations with other 
rare decays for a particular choice of some fundamental 
parameters (as in the case of small or large [S] tan [3 
in SUSY models with MFV). A precise determination of 
the decay rate of B^ would allow different NP 

scenarios to be disentangled. 

As shown in Table IH the present experimental limits 
on B^ are several orders of magnitude larger 

than SM expectations. Nevertheless, improved experi- 
mental bounds will restrict the allowed parameter space 
of several NP models. 

The search for the 5° t^t^ decay has been pre- 
sented in a previous paper 

In this paper, we present a search for B^ de- 
cay using data collected with the BaBAR detector fl3l| at 
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FIG. 1: Representative Feynman diagrams for ^ rr m 
the Standard Model. 



TABLE I: The expected branching fractions in the Standard 
Model 0] and the available upper limits (UL) at 90% C.L. 
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the PEP-II e+e" storage ring at SLAC. The collider is 
operated at the T(4S') resonance with asymmetric beam 
energies, producing a boost (/37 ~ 0.56) of the T(4S') 
along the collision axis. 

The dataset used consists of 384 x 10^ BB pairs 
accumulated at the T{AS) resonance ("on-resonance"), 
equivalent to an integrated luminosity of 347 fb"^, 
and 37 fb"^ accumulated at a center-of-mass (CM) en- 
ergy about 40MeV below the T(4S') resonance ("off- 
resonance"). The latter sample is used to characterize 
background contributions not originating from B decays. 

Hadronic two body decays of B mesons such as 
7r+7r~ and B^ — > K^n^ have the same event topology 
as the leptonic ones and are therefore the main source 
of background from B decays. We use Monte Carlo 
(MC) simulations [li] of B° ^ e+er, M+M", 
B'^ — > e^/i^ decays (signal) and B'^ tt+tt^ and 
B^ K^TT^ (background) of approximately 3 x 10^ 
events each to optimize event selection criteria and to 
estimate efficiencies. 

Charged particles are detected and their momenta 
measured by the combination of a silicon vertex tracker, 
consisting of five layers of double-sided silicon detectors, 
and a 40-layer central drift chamber, both operating in 
the 1.5-T magnetic field of a solenoid. The tracking 
system covers 92% of the solid angle in the CM frame. 
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Identification of charged hadrons is provided by the av- 
erage energy loss {dE/dx) in the tracking devices and 
by an internally-reflecting ring-imaging Cherenkov de- 
tector. For lepton identification, we also use the energy 
deposit in the electromagnetic calorimeter consisting of 
6580 CsI(Tl) crystals and the pattern of hits in resistive 
plate chambers (partially upgraded to limited streamer 
tubes for a subset of the data used in this analysis) inter- 
leaved with the passive material comprising the solenoid 
magnetic flux return. 

We reconstruct B'^ meson candidates from two oppo- 
sitely charged tracks originating from a common vertex. 
Signal events are characterized by two kinematic quanti- 
ties: 

mES = v'(s/2 + Po-Pb)V^^o'-p| (1) 

AE = E*B- (2) 

where \/s/2 is the beam energy in the CM frame, the sub- 
scripts and B refer to the initial T{AS) and to the B 
candidate in the laboratory frame, respectively, and the 
asterisk denotes the T(4S') rest frame. In Eq. the 
variable s is used as opposed to Eg because s is known 
with much greater precision and the resulting correlation 
between mEs and AE is nearly zero. For correctly recon- 
structed mesons, mEs peaks at the mass of the B'^ 
meson with RMS of about 2.5 MeV/c^, and AE peaks at 
zero with RMS of about 25 MeV. We require \ AE\ < 150 
MeV and toes > 5.2 GeV/c^. 

Since we use the pion mass hypothesis for the recon- 
struction of tracks, the distribution of AE peaks near 
zero for the tt+tt^ and leptonic modes and at -50 MeV 
for X±7r=F. The mass hypothesis does not affect the dis- 
tribution of Toes- 

Energy loss by electrons due to final state radiation or 
Bremsstrahlung in detector material leads to tails in the 
AE and toes distributions, in particular for the B'^ —> 
e+e" decay mode. We partially correct for this effect by 
adding the momentum of a photon emitted at a small 
angle from the track to the electron momentum. 

We apply stringent requirements on particle identifi- 
cation (PID) to reduce the contamination from misiden- 
tified hadrons and leptons. In this way, we retain 
~ 93% 73%) of the electrons (muons), with a mis- 
identification rate for pions of less than ~ 0.1% (~ 3%). 

According to the information provided by the PID, 
we separate our dataset into three samples, 2e, 2/i, and 
1/ile, containing events with two electrons, two muons 
and one muon and one electron, respectively. The rest 
of the dataset [h^h'^) comprises two oppositely charged 
hadrons and is used to characterize background contam- 
ination to the three signal samples. 

Based on MC simulations, we expect negligible cross 
feed of events between the leptonic and hadronic data 
samples. 



Contamination from non-resonant qq (q = u, d, s, c) 
and r+r~ production is reduced by exploiting their dif- 
ferent event topology with respect of that of the signal 
events. In particular, we examine the distribution of 
final-state particles in the rest frame of the T{AS) can- 
didate, in which the fragmentation of a BB pair (non- 
resonant event) produces an isotropic (jet-like) angular 
distribution of the particles. 

Non-Si? events ar e rej ected by requiring the cosine of 
the sphericity angle (l5!| to be |cos0s| < 0.8, and the 
second normalized Fox- Wolfram moment [lB| to be i?2 < 
0.95. In addition, we use a Fisher discriminant [l3| {^) in 
the maximum likelihood (ML) fit to separate the residual 
background from signal events. !F is constructed from 
the CM momentum pi and angle 9i of each particle i in 
the rest of the event (ROE) with respect to the thrust 
axis [18j of the B candidate. 

EE 0.5319- 0.6023^0 + 1.26987.2, (3) 

where Lq = Y^^^ Pi and L2 = Z^f^^ Pi ^os^ 6'i. The 
coefficients of the linear combination have been optimized 
on samples of signal and background simulated events. 
Since the variable T depends only on the ROE, we use 
the same coefficients for the three leptonic decays in the 
ML fit. 

The background from other BB events is found to be 
negligible after applying the PID requirements. Back- 
grounds originating from QED events (electrons and 
muons coming from e^e~ interactions) are rejected by 
requiring more than four charged tracks in the event. 

To ensure the quality of the measurement of the 
Cherenkov angle 9c, we require more than five detected 
Cherenkov photons and 9c > Q- For pion or lepton can- 
didates, in order to reject protons, we require 9c to be 
within 4(7 of the value expected for pions. For kaon can- 
didates, we require 9c to be within 4cr of the expected 
value for kaons. 

Applying the criteria described above, we select 67 
events in the 2e sample, 56 in the 2/i sample, 86 in the 
1/ile sample and w 94 x 10'^ in the h^h'^ sample. 

Among these events, the three signal yields are inde- 
pendently determined by ML fits on the 2e, 2/i and 1/ile 
samples. Each likelihood function is based on the vari- 
ables Toes 7 and . The probability density functions 
(PDFs) for the signal toes and AE distributions are pa- 
rameterized as: 

where /i is the maximum, CT/j/x, represent the standard 
deviation of the Gaussian component and cir/l describe 
the non-Gaussian tails of the PDF for a; > /i (R) and 
X < fi (L). The distribution for signal events is de- 
scribed by a Gaussian function with different RMS on 
the left and right side. The PDF of the background toes 
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distribution is parameterized by an ARGUS [l^ func- 
tion, the background Ai? distribution by a second order 
polynomial and the background T distribution by the 
sum of two Gaussian functions. Figure [5] shows the esti- 
mated background distributions for the three subsamples 
(solid lines) and, just for comparison, the corresponding 
signal PDFs obtained from Monte Carlo (dotted lines) 
with arbitrary normalization. 

We find that the residual background distributions of 
TOES; Ai? and T are the same in the three leptonic 
samples. This has been verified using data in the off- 
resonance sample and on-resonance events populating the 
kinematic sidebands (toes < 5.27 GeV/c^ or |Ai?| > 150 
MeV). 

In the fit the shape parameters for the — > l^V^ 
(signal) PDFs are obtained from the MC simulation with 
a correction factor that accounts for differences between 
data and MC, while the background PDF shape param- 
eters are determined on data with a procedure described 
below. 

We determine the parameters of the background PDFs 
by fitting their distribution on the h^'h'~ sample, where 
we use the Cherenkov angle to separate B'^ tt+tt" and 
SO ^ K^-K^. 

The yields of 5° tt+tt^ and K^t\^ in our h^h'^ 
sample are consistent with the results of the previous 
BaBAR analysis We find - 600 signal and - 3.5 x lO'' 
background events for tt+tt^, ^ 2200 signal and ^ 2.3 x 
lO'' background events for K^tt^. 

The background shape parameters in the i?" l^l'^ 
fit are fixed to the central values obtained in the fit to 
_B" 7r+7r~ and B^ K^w^ samples, and their er- 
rors are used to estimate the associated systematic un- 
certainty on the leptonic yields. 

We find no bias in the background shape parameters 
determined by the procedure described above on a large 
number of MC simulated h~^h'~ event samples. 

We correct for discrepancies between data and MC in 
the i?° — > /+/'^ signal shape parameters by rescaling 
the PDF parameters obtained from the simulation by 
the ratio between the values of the i?° vr+vr^ PDF 
parameters in data and MC. 

The knowledge of the rescaled shapes is limited by the 
size of the B'^ —> tt+tt" component in data, which causes 
a strong correlation among the parameters of each signal 
PDF. In order to avoid double-counting of these effects, 
we take the largest observed deviation as the systematic 
error induced on the leptonic yields. The errors on the 
signal yields due to the PDF shapes are ~ 1.1, ^ 0.4 and 
~ 0.2 events for the e^e^, /i^/i^ and e^/i^ channels, 
respectively. 

Our results are summarized in Table [ill We find no 
evidence of signal in any of the three modes. Using a 
bayesian approach, a 90% probability upper limit (UL) 



on the branching fraction (BF) is calculated as 

J C{BF)dBF / J £{BF)dBF = 0.9. (5) 
The BF is calculated as 



with Nw indicating the signal yield, ew the reconstruc- 
tion efhciency, and N^g the number of BB pairs in the 
dataset, N^g = (383.6 ± 4.2) x 10*^. We make the as- 
sumption that the T{AS) branching fractions to B~^B~ 
and B'^B^ are equal. 

The likelihood C{BF) is obtained by including in the 
likelihood function for the signal yield C{Nu>) the sys- 
tematic errors on Nui and the total number of BB pairs, 
and the statistical and systematic errors on the efhciency 
ew . We use the relation of Eq. ([6|) and assume Gaussian 
shapes for the errors. Figure [3] shows the likelihood dis- 
tributions of the three leptonic decays. 

We evaluate the efficiencies for individual selection cri- 
teria from MC simulation and correct the results for small 
differences between the simulation and the data. We take 
these observed differences as a measure of the systematic 
uncertainties on the efficiencies. 

The efficiency of PID requirements is calculated by us- 
ing MC simulations of signal events. It is then corrected 
with efficiency ratios computed on data and MC, as func- 
tion of track charge, momentum, and polar angle. We 
take into account the systematic error associated to this 
correction. 

The total systematic error on the efficiencies is ~ 4%, 
calculated as the sum in quadrature of all these contri- 
butions. 

In summary, we find no evidence of signal for i?" — > 
and place 90% confidence level upper limits on the 
branching fractions of i?^ — > e^e^, /i^/i^ and e^/i^. 

Table |TT] reports the efficiency, the number of signal 
events and the UL expected in each of these modes based 
on MC simulation for a sample of the size of our data 
sample. 

The present result on B^ e^/i^ and B'^ M^M^ 
improve the previous BaBAR upper limits |9] based on 
111 fb"^ 

The upper limit reported here for B'^ — s- e^e^ is higher 
than the value obtained in Q . In o ur p revious paper we 
used a largely frequentist approach [22] that does not ex- 
plicitly require a non-negative signal. The present results 
supercede our previous results: the analysis has a higher 
sensitivity, estimated from the value of the expected UL, 
and is based on a larger dataset that includes the sample 
used in 9]. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
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FIG. 2: Distribution of events in rriES (a,b,c), Ai? (d,e,f) and (g,h,i) for B —> e^e~ (left), B fJ-'^ IJ-~ (middle) and 
e*/i^ (right). The overlaid solid curves in each plot are the background sPlot [2^ distributions obtained when the 
corresponding component is ignored in the maximum likelihood fit and the likelihood maximized with respect to all the other 
components. The dotted line is the PDF obtained from signal Monte Carlo with arbitrary normalization. 



TABLE II: Efficiency (en')i number of signal events (Nui), 
90% C.L. Upper Limit on the BF {\JL{BF)) for the three lep- 
tonic decays 5° e^e~ , B^ — > iX^ and B^ — > e^p^ and 
the corresponding expected value based on MC simulation. 
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ew{%) Ng, \JL{BF) X 10-*^ Exp(UL) 

B°^e+e~ 16.6 ±0.3 0.6 ±2.1 11.3 7.4 

^ M^^" 15.7 ±0.2 -4.9 ± 1.4 5.2 5.9 

B°^e^/i^ 17.1 ±0.2 1.1 ±1.8 9^ 6.3 
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